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Abstract. As a consequence of continuous degeneracy of the order parameter of the superfluid 3 He quenched disorder in 
a form of aerogel gives rise both to a disruption of the orientational long-range order of the condensate and to a significant 
change of the order parameter itself. There exist a class of quasi-isotropic order parameters which are not sensitive to the 
disorienting effect of aerogel. While the BW order parameter belongs to this class the ABM does not. Possible candidate for 
the order parameter of the observed A-like phase is discussed. 



Average distance § fl between the strands of high poros- 
ity aerogels introduced as impurities in the superfluid 
3 He Lllyl does not differ much from the superfluid co- 
herence length ^o- This property is in a conflict with 
the condition of applicability of the conventional the- 
ory of superconducting alloys |3] n^Q 3> 1, where n is 
the average density of impurities. In aerogels n£? ~ 1 
and fluctuations of are significant. In a vicinity 
of the transition temperature T c effect of fluctuations 
can be described phenomenologically by introduction of 
random terms in the Ginzburg and Landau (GL) func- 
tional. For the scalar order parameter y/ such analysis 
has been made by Larkin and Ovchinnikov |4j]. Accord- 
ing to this analysis the most singular at T — > T c cor- 
rections to the order parameter originate from the term 
F n = N(0)Jri(r)\\i/\ 2 d 3 r (here N(0) is the density of 
states and r\ (r) - a random scalar), which describes local 
variations of T c . Generalization of this approach to the 
3x3 complex matrix order parameter An > is straightfor- 
ward. Additional term in the GL functional has a form: 



:N(0) 



Vji^AajA^r. 



(1) 



Now J]«(r) is a random tensor with zero average value. 
At the strength of t — > — t invariance tensor fj//(r) is real 

and symmetrical, its isotropic part Tjjj'(r) = \ r\mm{*)&ji 
describes local variations of T c — T c (r) due to fluctua- 
tions of the density of scatterers, it couples only to the 
overall amplitude of the order parameter A a /A 



Effect of T]jf (r) is similar to that for the scalar case 
- broadening of the transition and suppression of the 
absolute value of the order parameter. New properties 
originate from the anisotropic part T\u\*) = tj//( r ) — 
jrjn(r)8ji which describes local splitting of T c for differ- 



ent projections of angular momenta because of breaking 
of spherical symmetry by the aerogel strands. Contribu- 
tion of r}ji (r) to Fyf depends on the particular form of 
the order parameter and on its orientation with respect to 
the aerogel. Correspondingly rjj"\r) can both influence 
the form of order parameter and locally change its ori- 
entation. To make a qualitative guess of possible effect 

of Tlp\r) for different order parameters one can start 
with a uniform order parameter and substitute it into F^ . 
No essential difference with respect to the scalar case 

is expected if 7]^ A^jA*^ = 0. This is the case for the 

BW-order parameter of bulk B-phase: A™ = Ae'f R^j, 
where R a j - a real orthogonal matrix, but not for the 
bulk A-phase (ABM). Its order parameter has a form: 
a abm = (^/yfyd^fhj + ihj) with (Mj,hj) = 0, then 

r\f (r)A /li A* z ~ -r\fl]li £ 0, (here lj = e jlkmi n k ). This 
property of ABM-order parameter is of importance in 
connection with the identification of the A-like phase 
which has been observed in liquid 3 He in aerogel just 
below T c Hi]. 

Disorienting effect of the random tensor rjji (r) on the 
ABM-order parameter has been emphasized by Volovik 
Following argument of Larkin 1 8 ] and Imry and Ma 
he concluded that there may by no long-range or- 
der with the ABM-order parameter and suggested in- 
stead a glass state. In that state order parameter preserves 
ABM form locally but its orientation changes on a scale 
Lim ~ ^o/ 1 ? 2 - Characteristic value of rj 2 can be estimated 
e.g. from the measured width of the specific heat jump. 
This width is determined by the parameter a = r\ 2 / 
where T = (T — T ( )/T c . A region where fluctuations dom- 
inate is determined by the condition a ~ 1 |4[]. For 3 He 
using the data of Ref. lllOll for 22.5 bar one arrives at 



7] 2 ~ 1/5 and L /M ~ 5£ ~ 10" cm. 

Except for the disruption of the orientational order 

T]^ influences a form of the order parameter. For the 
scalar case only amplitude of the order parameter can 
change and its relative change is of the order of a <C 1 . 
For the triplet case the decrements of different compo- 
nents of A^j can be different. Moreover their change can 
be not small in case of degeneracy. For a triplet Cooper 
pairing there are altogether 18 extrema of the free en- 
ergy 1 1 1 ] and some of them can come close in energy. 
Contribution of fluctuations to the energy being small 
in comparison with the full condensation energy can be 
comparable or larger then the small difference of ener- 
gies of competing states. As a result fluctuations can shift 
the equilibrium and determine a choice of the order pa- 
rameter. Let us consider this possibility in more details. 
With account of the random term GL functional takes the 
form: 



rf 3 r[(T+iTj„(r))A w -A^- 



(T7;/(r) - - Wr)5 ;7 )A M; A* r 



dxj dxi 



K 3 



■ s=\ 
dxj dxi 



(2) 

where I s - 4-th order invariants in the expansion of the 
free energy over A^j. Coefficients Ai , ...As,^i ,K2,K$ - 
phenomenological constants. 

Because of a presence of r-dependent tensor TJ;/(r) 
in the functional (2) solutions of the corresponding GL- 
equations also depend on r. It is convenient to separate 
the order parameter into its average A^j and fluctuating 



part flii/(r) so that A„y(r) = A 



W + a W'( r )' Be low T c 
A^ij 7^ 0. The terms proportional to TJ«(r) in Eq. (2) 
dominate at T — > T c . Within a region |t| < T7 4 or a > 
1 a u,j(*) gives at least as important contribution to 
the physical quantities as A^j. It is difficult to analyze 
solutions of GL-equations within this region. Analysis is 
possible when a <C 1 and \a^j\ <C \A~uj\- In this region 
of parameters the iteration procedure of Ref.|4] renders 
equations for A^f. 



tK 



d 2 a 



dxfdxj 



-Animj, (4) 



Equations (3), (4) have the same structure as in the 
scalar case. Linear equation for a^j is solved by Fourier 
transformation and its solution has the following struc- 
ture: 



fl w -(k) 



{K 2 + k 2 ) 



jc~1/§. 



(5) 



Then for the average binary products entering Eq. (3) we 
have 



< T7 // (k)T7 m „(-k) > d 3 k 



(K 2 + k 2 ) 2 



(2k) 



(6) 

These terms in Eq. (3) give corrections of a relative order 
of <X <C 1 except for directions in which K = (Goldstone 
directions) when the integral diverges. This is the case 
for the increments of A^j and A* ■ at the infinitesimal 
rotation 9„: 



co^j = e n e jnr A^a;j = e n ei n % r 



(7) 



The divergency means that contribution of the aver- 
aged products of projections of a pi in directions deter- 
mined by Oifij to Eq. (3) is of a relative order of a^, 
X < 1. The value of X depends on a mechanism of 
cutting-off the divergency. The mechanism is provided 
by anharmonicity. Previously suggested evaluation of the 
diverging integral from the condition \api\ ~ \Api\ II 211 
overestimates it and renders X — 0. Better estimation is 
obtained by keeping anharmonic (third order) terms in 
Eq. (4). The obtained nonlinear equation can be treated 
within a mean-field approach, i.e. the third-order terms 
are substituted as aaa — > a < aa >, where a stays for any 
of components of a^j and only diverging binary averages 
< aa > s i„g are kept. Thus obtained linear equation has 
solution of a symbolic form (all indices are suppressed): 



A[(^ k) 2 + p<aa>]' 



(8) 



where A ~ l/T 2 is a typical value of /3 -coefficients in the 
GL-expansion and Q is a tensor Q r i = A^A*^ + A^/A* r . 
For binary averages in the Goldstone directions we have 



1 5 dl, 1 , d 3 I s 
N + 2 £ ft [ dAjj + 2 ( dA^jdAvndApi <^>' 



d 3 Is 



dA;jdA* vn dAp, 
and a^j(r): 



<a*vn a pi>)]=- < a f iiriij> ■ (3) 



< aa > s ing= const 



Q 2 f < 7777 > k 2 dk 
~&J [(^k) 2 + P<aa >~gf 



(9) 

This is an equation for < aa > S i n g- Its solution has a form 

Q 2 < 7777 > 



3/2 

< aa > sing = const. 



(10) 



1 5 d 2 I s d 2 I s 

la sl j + - £ As L - J. — a v „ + — ^, a* 



•W ■ 2tr i ^ } dA 



VII 



dA^dA^r- 1 



Here < Tjr/ >o= / < 77(0)77(1-) > d 3 r. With the aid 
of Eq.(8) < aa >«w can be estimated as < aa >j,w~ 



(AQ 2 a) 2 ^. At a — > the singular terms in the GL- 
equation give a principal part of the contribution of fluc- 
tuations to this equation. GL-equation now has a struc- 
ture: 



8A 



-[F 



(0) 



+F fi. 

1 reg 



F fl 1 = 0. 



(ii) 



Consider a set of competing states which are nearly 
degenerate with ABM: 



Fc-F, 



ABM 



Fabm 



= e < 1, 



(12) 



then, in Eq. (11) the first term is of the order of e, the 
second - of the order of a and the third - of the order of 
(a) 2 / 3 . A character of solution depends on a relative val- 
ues of small parameters e and (a) 2 / 3 . If e 3> (a) 2 / 3 the 
order parameter is determined as a minimum of the un- 
perturbed GL-functional (supposedly ABM-state) 
and fluctuations give small corrections to its form. In the 
opposite limit e <c (a) 2 / 3 the form of the order parame- 
ter is determined by the fluctuations: 



SA 



-ft 



f< ] 



0. 



This equation is satisfied if Q 



ji 



ist.8:, 



or, equivalently: 



,«<> 



0. 



;i3) 



(14) 



(15) 



As has been discussed above, for the order parameters, 
meeting condition (15) there is no disorienting effect 

of the random tensor TJ an d long-range order is pre- 
served. Generally this condition specifies a degenerate 
family of so-called "robust", or quasi-isotropic states. For 
partial lifting of the degeneracy the remaining terms in 
the free energy have to be minimized 



8A 



-[f^+fB=o 



(16) 



with Eq. (14) as a constraint. Eq.(16) determines in par- 
ticular the overall amplitude A. It has to be verified then 
that the found solutions satisfy condition (12). If it is 
the case there exist a temperature region where condition 
£ -C (ct) 2 / 3 is met as well since £ does not depend on 
temperature but a grows as l/y^frf at T — > T c . Within 
this temperature region the above procedure of finding 
solution of GL-equation which are determined by fluctu- 
ations is justified. 

In case of A-like phase the competing states have to 
be of equal spin pairing (ESP) type. It follows from the 



experimental observation that magnetic susceptibility of 
the emerging phase coincides with that of normal and 
ABM-phases 1 5 ] . A general form of the order parameter 
for ESP-state is: 



,ESP 



: ^Tfi^imj + inrf + e^lj + ipj)], (17) 



where and are mutually orthogonal unit vectors, 
nij , tij , lj , pj - arbitrary real vectors. The order parameter 
A^f meets criterion (14) if vectors mj,njjj,pj satisfy 
the equation 



nijini +tijtii + l jk + pjpi = Sji . 



(18) 



A detailed description of this solution is presented else- 
where 1 13]. If additional constraint of the absence of 
spontaneous magnetization is imposed the obtained so- 
lution can be rewritten as: 



-icosycy) + + isinycj)], 

(19) 

where dj,bj,cj - mutually orthogonal unit vectors, y/ 
and y- parameters. 

Now the condition (12) has to be verified. For the order 
parameter Eq. (19) e is expressed in terms of coefficients 
Pi as: 

e= ft3-4fa (20 ) 

9^2 + ^3+5^45' 1 j 

where /3n = /3i + ft etc. For the weak coupling values 
of j3 -coefficients e = 1/19 i.e. small, but in that case 
BW-phase is more favorable then ABM. The most im- 
portant are strong coupling corrections to the combina- 
tion /J45 since the weak coupling value ^45 = deter- 
mines a boundary between ABM and axiplanar phases 
in the space of j3 -parameters 1 14]. For the moment it is 
difficult to make definite conclusions about the value of 
[$45 realized in liquid 3 He. Combinations of ft , ...ft de- 
duced from the analysis of thermodynamic data I il5ll do 
not form a complete set and can not define e without am- 
biguity. On the other hand the values of ft , ...ft found 
theoretically from model assumptions 11611 do not agree 
with the experimental findings in particular about ft. 

Eq. (19) gives solution of GL-equation of a zeroth or- 
der on the small parameter e/a 2 / 3 . In the real 3 He this 
ratio can be not very small so that the finite deviation 
from "robust" state occur and disorienting effect of the 
random tensor rjji (r) comes into play. The outlined pro- 
cedure offers in that case only a local form of the order 
parameter. Principal contribution to the integral Eq.(9) 
comes from the wave vectors (^o^) 2 ~ P < aa > S ing, or 
length scales L c -o// ~ % /(ccq 2 ) 1 ^, here q 2 — Q 2 /A 2 is 
the normalized deviation of the order parameter from 
the "robust" form. The local form of the order param- 
eter is preserved if it changes its orientation in space 



on a length scale which is much greater then L c _ yy. 
This is the case for Imry and Ma disorienting effect: 
LiM/L c - ff ~ 1 /(ag 2 ) 2 / 3 » 1. The resulting glass-type 
state is based on a "nearly robust" state and the charac- 
teristic size of the "domains" is much greater then in the 
case of the ABM-based glass state. The most important 
difference between the two states is that the order pa- 
rameter A^j unlike the A^ B j M does not have the combined 
gauge-orbital symmetry. As a consequence the property 
of superfluidity for A^ ; - is preserved even in the glass- 
type state. Formally it means that the averages of a type 
< A^ / (r)A^ / (r) > are finite and can be used as the order 
parameter of the "superfiuid glass" state 11711 . 

So, one can conclude that at a triplet Cooper pairing 
quenched disorder of a type of "fluctuations of the tran- 
sition temperature" rjji (r) influences both form of the or- 
der parameter and its local orientation. There are "robust" 
classes of order parameters, for which disorienting effect 
disappears. If the form of the order parameter required 
by the minimization of the initial free energy is not "ro- 
bust" there is a mechanism of adjustment which tends 
to bring the order parameter closer to the "robust" form 
and to decrease the disorienting effect of the random ten- 
sor r\ji (r). Both disorienting effect and the mechanism of 
adjustment originate from degeneracy of the unperturbed 
free energy with respect to orbital rotations of the order 
parameter. 
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